An ion ring in a linear multipole trap for optical frequency metrology 
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A ring crystal of ions trapped in a linear multipole trap is studied as a basis for an optical 
frequency standard. The equilibrium conditions and cooling possibilities are discussed through an 
analytical model and molecular dynamics simulations. A configuration which reduces the frequency 
' sensitivity to the fluctuations of the number of trapped ions is proposed. The systematic shifts for 

04 I the electric quadrupole transition of calcium ions are evaluated for this ring configuration. This 

study shows that a ring of 10 or 20 ions allows to reach a short term stability better than for a 
' single ion without introducing limiting long term fiuctuations. 
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PACS numbers: 37.10.Ty (Ion trapping) 37.10.Rs (Ion cooling) 06.30.Ft (Time and frequency ) 



^ ■ I. INTRODUCTION 



Progress in atom and ion cooling and trapping, laser stabilization and high-resolution spectroscopy makes narrow 
optical transitions the best candidates for frequency standards. Three experimental systems take the optical clocks 
^ ■ to their highest performances: a laser-cooled neutral atom cloud in a MOT [l|, Q, an ensemble of laser-cooled neutral 
Q^l atoms in an optical lattice and single trapped laser-cooled ions [6- 9] . All systems have advantages and drawbacks. 
The main advantage of single ion experiments is the quasi-perfect control of the internal and external degrees of 
^ freedom, permitted by the trapping in radiofrequency (rf) electric fields in the quadrupole configuration. This ultimate 
^ . control, possible only with a single particle, is paid for by a low signal-to-noise ratio which requires a long integration 
' time, resulting in a small short term stability when compared to neutral atom optical clocks where 10^ to 10* atoms 
\^ , are interrogated at the same time. 

In the same time, microwave frequency standards with large clouds of trapped ions (> 10*^ ions) [l3[ll| are developed 
[ to assure a good short term stability (characterized by an Allan deviation a-y{T) < 10~^^ / y^) and an excellent long 
^f) • term stability cry{T) < 10~^^ per day). In the context of deep space navigation [12], precision is not a major issue 
' and a fractional frequency uncertainty of 10~^^ is admitted to be sufficient. The main systematic effect which limits 
. the long term stability is the second order Doppler effect induced by the rf driven motion (called micromotion). This 
' effect depends on the size of the cloud through the maximum amplitude of the electric field seen by the ions and 
^ , requires a high stability of the ion number over a long time scale. For an equal number of trapped ions, this shift is 
reduced by using linear multipole traps where the electric field amplitude is almost null in the center part of the trap 
[isj . Large ion samples can thus be trapped with reduced micromotion, compared to the quadrupole geometry. For 
• a dense enough ion cloud, one can show [12] that the maximum second order Doppler shift SfD2 = — fo (i'^)/(2c^) 
simply writes 

0JD2 — -JO-X 77 TT-^, (Ij 

87reom(/c — l)c^ 

where /o is the frequency of the atomic clock transition, is the number of ions per unit length, m is the mass of 
the ions, q their charge, 2fc is the number of electrodes in the multipole trap and c the speed of light. The (fc — 1)^^ 
scaling justifies the choice of a higher order multipole [T^ for challenging clock applications. Indeed, in [l2| a 16-pole 
trap is used to trap ions during the interrogation of their clock transition. Multipole rf traps are also widely used to 
study cold reactive collisions [15] as cooling with a buffer gas from 300 K to 4 K 16] allows to control the kinetic 
energy of the collisions and to reproduce astrophysical conditions. 
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Laser cooling and crystal structures of cold ions in multipole traps are far less studied than in quadrupole traps. 
Laser cooling and observation of ion crystals in a linear octupole trap have been reported in (TtI [18| . A semi-analytical 
and numerical study of the structure, scaling laws and phase transitions of cold ions in an isotropic 3D octupole can 
be found in These two studies show that a cold ion cloud can be described as a hollow core system, resulting 
from the balance between the Coulomb repulsion and a trapping potential nearly flat in the center and very steep at 
the border. Such a geometry can be also deduced by a cold fluid model [lO]. When the number of trapped ions is 
reduced to the order of 100 or less, simulations show that, for certain trapping parameters, the tube formed by the 
ions reduces to a ring crystal centered on the symmetry axis. 

In this article, we propose to use laser-cooled ions trapped in a linear multipole rf trap and organized in a ring 
structure as a basis for an optical clock. The aim is not to compete with the highest precision a single ion optical clock 
can offer, but to propose a trade-off where a somewhat lower precision is compensated by a gain in short term stability 
offered by the interrogation of several ions at the same time. The main technological challenge in the realization of 
a single ion frequency standard is the clock laser which has to reach frequency stabilities of the order of 1 Hz/s. 
The present proposal allows to reach total clock performances comparable to a single-ion experiment by relaxing the 
constraints on the laser performance by at least one order of magnitude. The improvement can be illustrated by the 
Allan variance which is used to quantify the stability : 



with Q the quality factor of the measured transition, (SNR) the signal-to-noise ratio of the excitation probability 
measurement and Tc the cycle time required to build an atomic signal to counteract on the laser frequency. The 
(SNR) is limited by the quantum projection noise [2T| to \/N and using TV = 10 or 20 ions allows to reach the same 
stability value 10 or 20 times faster than a single ion clock, if the line is not broadened. 

Compared to a chain of ions in a linear quadrupole trap, this ring configuration has the advantage of having the 
same symmetry as a transverse laser beam intensity profile assuring that every ion sees the same laser intensity. The 
second advantage of this symmetry is that the motion along the trap axis is characterized by a single oscillation 
frequency. Furthermore, the trapping parameters can be chosen to constrain the radial size of the ring, independent 
from the number of trapped ions, to first order (if this number obeys some stability conditions). As a consequence, the 
possible loss of one ion during long time operation induces only a second order perturbation on the ring equilibrium 
radius and causes very small frequency fiuctuations. The compactness and the symmetry of a ring compared to a 
chain made of the same number of ions is another advantage regarding the dispersion of the systematic shifts induced 
by local electric or magnetic fields. 

The planar crystal structure of the present study could also be of use for quantum computation experiments. Planar 
crystals in a quadrupole trap are proposed in 22] as an appropriate system for large scale quantum computation and 
their structure is studied in (23j by molecular dynamics simulations. The ring configuration we propose here for 
trapped ions could also be used for quantum simulations like demonstrated in [24] with several Rydberg atoms at 
each site, to investigate many-body quantum states. 

In the following section, we start by describing the features concerning the trapping and cooling of the sample, 
relevant for an optical clock of this kind. We then identify two configurations that obey the clock operation conditions: 
a ring of 10 ions and one of 20 ions. In section Hill the scaling laws and numerical values of the systematic shifts due 
to the ring configuration are given. In a second step (section ITV)) . the systematic shifts usually encountered by a single 
ion clock are estimated for this ring configuration. In the last part, all contributions are summarized to propose an 
uncertainty budget for the 10- and for the 20- ions ring clock. The numerical calculations quantifying all the effects 
use the optical quadrupole transition of '^^Ca^ as the clock transition [251, 126^] and are carried out for a linear octupole 
trap (2k = 8). However, most of the equations can be transposed to other species and to higher order traps. 



In the adiabatic approximation, the static pseudopotential associated to the rf electric field allows to explain the 
dynamics and equilibrium positions of the ions in the transverse plane of the trap. For a perfect linear multipole trap 




(2) 



II. SELF-ORGANIZATION OF IONS IN A RING 



of order 2k, it is defined like 




(3) 



where r is the ion distance to the trap axis, Vq is the amplitude of the rf potential difference between two neighbouring 
electrodes, r2/27r the frequency of this applied voltage and tq the inner radius of the trap. To effectively trap the 
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FIG. 1: Solid line: effective static potential (in eV) in the radial plane of an octupole trap {2k — 8), resulting from 
the contributions of V*{r), Vstat- Dashed line: the harmonic approximation of the potential well defined by r,„i„ 
(Eq.O and We// (Eq.[5]). The potentials are defined by u}:,/2'k ^ 1 MHz, nj2n = 20 MHz, Vb = 394.4 V and 

ro = 200 ^m, for calcium ions (to — 40 a.m.u). 



ions, a static potential has to be applied on electrodes placed on both ends of the linear trap. This geometry creates 
a potential well along the axis of symmetry Oz which we can consider as harmonic in the centre of the trap. The full 
contribution of this static potential can be written like 



1 / 

^ .2 ( „2 ' 



Vstat = ^m^i [' 2) 

where the harmonic well is characterized by the oscillation frequency w^. When the multipole is not a quadrupole 
{2k 7^ 4), the two contributions V*{r) and Vstat for the trapping potential result in a potential shape (solid line in 
figure [T|) , where the minimum is shifted from r = to r = rmin defined by 



2fc-4 _ 1 f 2mnu!zrQ 



1 V kqVo 



(5) 



A. Equilibrium 

The equilibrium position of a set of ions results from the balance between the Coulomb repulsion and the trapping 
potential. An analysis of the forces along the radial direction shows that two regimes of organization can be defined, 
based on the strength of the deconfinement contribution — mw^r^/4 compared to the Coulomb repulsion induced on 
an ion at radial distance R by all the other ions. In the assumption of a ring of radius R in a single plane, the Coulomb 
potential energy for a set of N ions Ep'{N) is 

E-iN)^^^^ (6) 
P ^ ' ineo 2 2R ^ ' 

with Si{N) = J2n=i 1/ sin(7m/A'^). If this contribution is negligible compared to the deconfinement, the radius of 
equilibrium R does not depend on the number of ions TV. In the prospect of a clock where the major systematic shifts 
depend on the equilibrium position R through the local rf electric field, this condition has to be fulfilled to allow the 
ion number not to be strictly reproduced from day to day. In the following, we consider that the configuration is such 
that, at first order of approximation, the equilibrium position R is given by the minimum potential rmin defined by 
Eq.[5]and does not depend on the number of ions in the ring. The Coulomb repulsion can be treated as a perturbation 
and the radial position shift e that it induces depends on various trap parameters like 

^ (7) 

4^60 4(fc - 2)r2„„mw2 
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With the ion number N and ring radius considered in the following, the relative Coulomb shift e/r„ 
0.02, justifying the assumption that R ~ r^irt to deduce the scaling laws for the systematic shifts. 

With the typical trapping parameters used through this article, local potential depths equivalent to 10 K can be 
made while temperatures as low as 10 mK can be reached in the radial direction by Doppler laser cooling. Laser 
cooled ions will settle in the bottom of this well and at this level, the radial pseudopotcntial can be approximated by 
a harmonic potential centered on rmin (see Figure [l}. An analytical analysis shows that the second order expansion 
of the radial potential around rrmn, V*{r) + Vstat{r) ~ V*{rmin) + Vstat{rmin) + muj^j: j:{r - rmin)'^ /2, is characterized 
by an effective oscillation frequency Weff which depends only on the axial frequency and the order of the muitipoie 
by 



= Vfc - 2 LJ^. (8) 



This shows that the strength of the axial trapping ujz also defines the strength of the radial trapping around the 
equilibrium position R. 

To make relevant predictions based on the pseudopotential (Eq|3]), the adiabatic approximation at position Tmin has 
to be satisfied. Contrary to quadrupole traps, there is no absolute criterion for adiabatic operation of the muitipoie 
trap. However, the Mathieu parameter used in quadrupole traps to characterize the trajectories of ions can be 
gen eralized to muitipoie with the main difference that this new parameter rj depends on the ion's location in the trap 
[l^H^I. For a perfect muitipoie, this parameter is defined like 

vir) = Hk~l)-^(^y ' (9) 



mQ'^rQ \rQ 



and has been empirically limited to 0.3 to guarantee stability of the trajectories 15|. A more recent experimental 
study of the loss mechanism in a 22-pole trap [1^ complementary to a model of effective trapping volume [l^ 
has demonstrated stability up to ri{r) < 0.36 ± 0.02. Moreover, clock operation requires also a small micromotion 
amplitude to limit all the systematic effects induced by this motion and to minimize the rf heating that may occur, 
as in quadrupole traps. For this reason, and by analogy with linear quadrupole traps, we limit the set of trapping 
parameters to keep rj < 0.2. For a ring of ions at a distance R = Tmin from the center and including Eq. [5]into Eq. |9l 
the adiabatic parameter takes the simple form 



ry(i?) = 2^/fc^^ (10) 

which does not explicitly depend on the equilibrium position nor on the number of ions, and immediately fixes the 
range of the trapping frequency V,/2tt once the axial oscillation frequency is chosen. Indeed, in an octupole trap with 
uJz/2tt = 1 MHz, it takes ri/27r > 17 MHz to assure ij < 0.2. Furthermore, the amplitude of the micromotion 
also scales like ri{R) but increases with the distance to the center like 



1 

'2(fc~ 1)' 



5R^^Rj^{R)— — ^. (11) 



Therefore, a radial equilibrium position as small as possible is chosen, to reduce the rf induced heating as well as 
the rf induced systematic shifts, which depend on the local rf electric field and scale like R^~^. The lower limit of 
this radius is set by the ion-ion distance. For a distance too small, the one ring configuration is not stable and the 
system relaxes to a double ring where the sum of the Coulomb and trapping potential energies is lower. Simple energy 
minimization, confirmed by molecular dynamics (MD) simulations, shows that in a double ring of an even number of 
ions, the ions alternate from one ring to the other in a configuration analog to the zig-zag chain already observed in 
a quadrupole linear trap (30| (see Figure [5]). Writing the equilibrium condition along Oz for a double ring of N ions 
requires that R should be smaller than a limit Ri, which sets the lower stability limit for a one-ring configuration. 
This limit size depends only on the strength of the axial potential and the number of ions as 



1/3 



1/3 



and therefore scales like N. 
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FIG. 2: (Color online) Equilibrium position of 20 calcium ions around the trap symmetry axis, calculated by MD 
simulations in the pseudo-potential of a linear octupole trap. The potential is defined by uiz/lir— 1 MHz and il^/27r 
= 20 MHz and for 20 ions, the stability limit is Ri = 23/im. Left: Vo ~ 3142 V and r^— 400 /xm, resulting in 
f-min > Rh the One ring configuration is stable (R — 28 /im). Right: Vq ~ 5771 V and — 400/im, resulting in 
'^min < Rii the stable configuration is a double ring analog to a closed zig-zag chain. Here the separation between 
the two planes is of the order of 3 im\ and their radius is 21 /im. 



B. Doppler cooling 

So far, we have been concerned by the static properties of the ring configuration. In this second section, we 
focus on the thermodynamic issues related to ion temperature required for clock operation. We first introduce these 
requirements and their consequences for the geometry of the set-up and then propose two configurations that fulfill 
these requirements. 

When the trapping parameters are adjusted to obtain a single ring of ions lying in the z = plane, it seems 
mandatory to propagate the clock laser along the Oz direction, such that no line broadening is introduced by dispersion 
of the excitation strength on the clock transition from one ion to another. In this configuration, there is no micro- 
motion along the laser axis and the first order Doppler effect is sensitive only to the thermal motion along the axial 
direction (called macro-motion). To compete with the existing clocks, this effect has to be canceled. The method used 
for trapped ion frequency standards is to set the ion motion in the Lamb-Dicke regime where the phase modulation 
induced by the ion motion inside the clock laser wave is small compared to 27r. For optical transitions, the wavelength 
is so small that laser cooling is required to reach that stage. One of the major issues of laser cooled ion trap design 
is to make sure that the Lamb-Dicke regime can be reached by Doppler laser cooling so that no other laser cooling 
process is required for clock operation [31| . 

When the ion motion can be described by an harmonic oscillation Zcoswzi, a classical expansion of the laser-atom 
interaction (32j shows that the Doppler spectrum can be understood like the superposition of sidebands at frequency 
/o ± nLOz/^TT with amplitude proportional to | J^(fcLZ)|, if the transitions are not saturated. J„ is the Bessel function 
of order n and the modulation index k^Z depends on the clock laser wavevector and the amplitude of oscillation 
Z = Vz/uJz- The velocity amplitude of the oscillation Vz depends on the temperature reached by the cooling Tj. 
When the Doppler laser cooling limit is reached, Vz = \J jm where 7 is the spontaneous decay rate on the cooling 
transition (at 397 nm for calcium ion). For the Ca^ clock transition at 729 nm with ions at the Doppler limit Tj^i^ 
Td = 0.54 mK, with I/7 = 7 x 10~^ s), the modulation index is 27r x 6.3 x 10^ /uj^- Choosing tj^/27r = 1 MHz is 
sufficient to assure a modulation index smaller than 1 which sets the motion in the Lamb-Dicke regime and leads 
to a central band nearly 10 times stronger than the first oscillation sidebands, which is a good enough operation 
condition for an optical clock[2^. We will see later in the text that further increase of the axial oscillation frequency 
is detrimental to the performances of the clock and the chosen value appears to be a good compromise between a 
small oscillation amplitude and small systematic shifts. 

All the considerations introduced in the previous part about the equilibrium position of the ions inside the pseu- 
dopotential are not sufficient to control the dynamics of the set of ions. To get more insight in the system, molecular 
dynamics (MD) is used to simulate the ion motion under Coulomb repulsion and laser cooling, in the full rf potential 

$;,(r, t) = Vq/2 cos(m)(r/ro)'= cos(/ca) (13) 



where (r, a) are the polar coordinates in the (a;, y) plane. Details of these simulations will be published elsewhere 
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[33| . These simulations show that in the single ring configuration, the axial and radial degrees of freedom are very 
weakly coupled and the thermal equilibrium of each degree is characterized by very different temperatures. With 
laser cooling on the three directions of space, the simulations show that for a decoupled system, the Doppler limit 
temperature is reached in the axial direction and the temperatures in the radial plane {Tr,Ta) reach 10 mK in the 
best configuration (laser detuning set at —7/2 and coupling strength defined by a Rabi frequency equal to 7/2). The 
MD simulations show that if the trapping parameters result in a radius of the ring too close to the stability limit 
Ri, the axial and radial degrees of freedom are coupled and the temperature for the axial motion does not reach the 
Doppler limit. Reaching the Doppler limit in the axial direction is mandatory for the Lamb-Dicke regime, and we can 
design potentials able to trap a ring made of 10 to 100 ions in a single plane, fulfilling this condition. 

For clock operation additional conditions have to be fulfilled. Indeed, to prevent light-shift on the clock transition, 
the cooling lasers must not be applied onto the ions during the clock transition excitation . During this excitation, 
which needs to last at least several milliseconds, the system must remain at the same temperature. With simulations, 
we have found trapping and cooling conditions compatible with these requirements only for rings smaller than 20 
ions. The 100-ions ring heats up very fast and, in the best configuration, the axial temperature increases from 0.54 
mK (the Doppler limit) to 2 mK in 1 ms. This rapid increase is due to rf heating as the micromotion amplitude 
increases with R (see Eq. pTjl ) which increases like N. On the contrary, for 10 ions set in a ring of radius 20 /xm 
or for 20 ions set in a ring of radius 40 /im, the cooling laser can be switched off for several ms with very little 
heating. To reach good statistics we build a sequence of 2 ms of dark time (the clock interrogation) followed by 2 ms 
of laser cooling, and measure the temperature throughout the dark times. For the 10-ions ring, after 120 periods of 
cooling/non cooling time, the average axial temperature is 0.47 mK and the FWHM of the distribution is 0.32 mK. 
Clock operation conditions for a 10-ions ring can thus be fulfilled (similar results can be obtained for the 20 ion ring). 
For comparison, all the configurations have the same axial and rf frequency u!z/2tt = 1 MHz, f2/27r = 20 MHz, and 
therefore the same adiabatic coefficient at the position of the ions ?y(i?) = 0.17. The ring radius was controlled by the 
amplitude of the rf electric field. 

The constraints on the trapping parameters and their relation can be summarized as follows: the minimum value 
of ujz is fixed by the Lamb-Dicke regime. Once ujz chosen, the value of f2 is fixed by the adiabacity criteria (Eq. llOp . 
Then, for a given number of trapped ions N, the ratio Vq/vq is chosen to match r„iin with the smallest R that allows 
to keep the ions in a single ring (i? > Ri) which does not heat up when the cooling lasers are off. The MD analysis 
shows that it is possible to form a stable ring of 10 to 20 ions, with a temperature close to the Doppler limit in the 
axial direction and of the order of 10 mK in the radial direction. In the second part of this paper we evaluate frequency 
shifts on the clock transition for two configurations tested by MD simulations : a 10-ions ring of radius 20 fim and 
a 20-ions ring of radius 40 /xm. With the chosen trapping frequencies, it takes Vq = 394 V for ro = 200 /xm to reach 
R = 20/im and Vq = 1578 V for tq — 400 /xm to reach R = 40/xm. Furthermore, we assume an axial temperature of 
0.54 mK (the Doppler limit Tjj for calcium ion) and a radial temperature of 10 mK. This temperature controls 
the thermal motion amplitude along the radial direction SR through 



and is responsible for line broadening through every systematic shift depending on the radial position R. To give an 
order of magnitude of this effect, for calcium ions in a linear octupole trap with ujz/'2n=l MHz, the thermal motion 
amplitude in the radial direction SR is 0.23 /xm, to be compared to the amplitude of the micro-motion which 
is 0.56 /im for a 20 fim radius (and 1.1 /xm for a 40 fim radius). These equilibrium and dynamic properties control 
several systematic shifts which are discussed in the following section. 



Compared to a single ion experiment, additional systematic shifts are introduced by the fact that the rf electric 
field is not null at the ion's equilibrium position. In this section, we focus on these extra effects and highlight their 
scaling law to point out the best compromise for a clock operation. The performances of a clock are evaluated by its 
accuracy and stability and the laws scaling their dependance with the number of ions are different. For every shift 
depending on the equilibrium position R (which scales with TV to obey the equilibrium condition R > Ri), there is 
a line broadening induced by the thermal oscillation amplitude SR (which does not depend on TV) and a long term 
fluctuation induced by the possible loss of one ion. This last effect depends on Ae = e(TV) — e{N — 1) which scales 
with l/N"^. All shifts are evaluated for the quadrupole transition at 729 nm (Si/2,T\^j — ^h/2i^j) of a ""^Ca^ ion 
which has no hyperfine structure. The trapping frequency is chosen to be ^1/2tt — 20 MHz, with an axial confinement 
of wj27r = 1 MHz. 




(14) 



III. SYSTEMATIC SHIFTS INDUCED BY THE RF ELECTRIC FIELD 
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A. Micromotion induced second order Doppler effect 

The second order Doppler shift induced by the micro-motion is known to be the hmiting factor for microwave ion 
clocks in multipole traps. The velocity amplitude of the micromotion is equal to flSR^^ and the shift induced on 
the clock transition frequency /o, can be expressed like 

In analogy to large clouds in a multipole trap (Eq. ([T])), this shift scales with l/(fc — 1), pleading for higher order 
multipole traps. But contrary to large clouds, the shift does not depend explicitly on the number of trapped ions but 
only on the trapping parameters, which can be very well controlled. The compromise to make for a clock operation 
is clearly visible from Eq. : the reduction of the second order Doppler shift asks for uJz and R values as small as 
possible whereas the ion ring stability and cooling set a lower limit to these two parameters. 

For a ring of radius 20/im (respectively 40/im) in an octupole trap {k = 4), the shift is 1.46 x 10~^* x fo = 6.0 Hz 
(respectively 5.85 x 10^^"* x /q = 24.1 Hz) on the 729 nm optical clock transition of Ca"*". The shift itself does not 
reduce the precision and stability of the clock, but its uncertainty and fluctuations can. 

The slow (on the rf time scale) thermal oscillation of the ions in the radial plane is responsible for a dispersion in the 
radial position ±dR and implies a broadening of the second order Doppler shift Adf]j2 — ±<5/D22(5i?/-R (full width), 
which scales like uJzR, which itself is proportional to N and is lower than one Hertz with the typical parameters we 
have chosen {ASfD2 = ±3.4 x 10'^^ x /o = ±0.14 Hz for R = 20iJ,m and ±6.7 x IQ-^^ x fo = ±0.28 Hz for R = 40nm). 

In long term operation, an ion can be lost and the shift in the radial equilibrium position Ae induces a modification 
of the frequency shift equal to SfD2'2Ae/R, which scales like 1/A^. For the 10 ion ring of radius 20/im, the fractional 
fluctuation is 8.5 x 10~^^ x /o = 0.03 Hz and it is 5.1 x 10^^^ x fo for a 20 ion ring of radius 40/im. One can see that 
the 10""'^^ threshold we aim at for the long term stability can already be reached with a ring configuration of 10 ions. 
Increasing the number of ions decreases the long term fluctuations. 

B. Stark effect induced by the electric trapping fields 

The Stark effect results from the contribution from a scalar shift, which is proportional to the square amplitude of 
the electric field -E'^/Z^ ± ^dc (^^^ ^/-^ being due to the time averaging of the rf component) and a contribution from 
a tensorial part, which is sensitive to the angle between the magnetic field and the electric field at the position of 
the ion [sjj. The scalar and tensorial shifts are proportional to the differential polarizability between the upper and 
lower states of the clock transition and the scalar contribution Aa^ = —1.1 x 10~^ Hz/(V/m)^ [s^ is independent 
of their Zeeman sublevels. For calcium ions, the tensorial contribution is of the same order of magnitude than the 
scalar one [1^ [H, [s^ and depends on the Zeeman sublevels Mj chosen for the D5/2 state (the S1/2 state has no 
contribution to the tensorial shift): Aa^ = -6.1 x lO"'^ x f{Mj) Hz/(V/m)2 [35I with f{Mj = ±1/2) = -4/5, 
f{Mj = ±3/2) = -1/5 and f{M,j = ±5/2) = 1. 

The amplitude of the rf electric field is deduced from Eq. (jl3p and for ions at distance R from the center, matching 
the potential minimum rmim its contribution to the scalar Stark shift can be written like 

The scalar Stark shift induced by the dc electric field can be expressed with respect to the rf contribution : 

SfUdc) = 5/0(r/)i^^^ = Sflirf)^. (17) 

This relation clearly shows that for ri{R) < 0.2, the contribution from the dc field is negligible compared to the rf 
field contribution. As the Coulomb repulsion is small compared to the dc radial electric field, the Stark effect induced 
by the Coulomb interaction is negligible compared to the dc contribution, which is itself negligible compared to the 
rf contribution. Therefore, only the rf contribution is taken into consideration in the uncertainty budget. 

Like the second order Doppler shift (Eq. [T5|) . the Stark effect scales with l/(fc — 1) and is proportional to uj^R^. 
The new parameters that are introduced in the shift dependence are the mass of the ion and the trapping frequency. 
With the same parameters as mentioned before, Sfg{rf) = 1.0 x 10""'^^ x /o = 4.1 Hz and the line broadening induced 
is ±2.3 X 10^^^ X /q ±0.09 Hz for the ring radius of 20 fj,m. For the ring radius of 40 /im, the shift Sfg{rf) is 
4.0 X 10^^* X /o = 16.5 Hz and the induced broadening is ±4.6 x 10~^^ x fo = ±0.19 Hz. The fractional long term 
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instability induced by AN = —1 is 5.8 x lO"""^^ for the 10 ion ring of radius 20 /im (and 3.5 x 10""'^'' for the 20 ion 
ring of radius 40 /xm) . 

The tensorial Stark shift induced by the rf electric field depends on the angle 9 between the local electric field and 
the magnetic field like [13] 

Sflirf)^~lAayiM,f-^^2fl^^, (18) 

In an rf trap, the amplitude of the rf electric field is constant over a ring but its direction rotates in the radial 
plane. The projection of this electric field on a static frame is defined by (cos(fc — l)a, sin(fc — l)a) if a is the angle 
between the ion radial direction and a reference axis chosen in the radial plane. If the geometric configuration is 
such that 6 varies from one ion to the other, the dispersion of the Stark shift induces a line broadening of the order 
of 1 Hz for the 20 /xm radius ring and of 4 Hz for the 40 /im radius ring. Such a broadening, which is bigger than 
any other broadening calculated so far, can be avoided by choosing a magnetic field oriented along the trap axis 
(6* = 7r/2), to assure the same Stark shift for all the ions of the ring. In this configuration and with the same trapping 
parameters, the tensorial Stark effect is Sfg{rf)/f{Mj) = —0.3 x 10^^'' x /o = —1.1 H, with a line broadening 
±0.6 X 10-16 X fo = ±0.02 Hz, for the ring radius of 20 /im {6ps{rf )/ f{Mj) = -1.1 x 10"" x /o = -4.6 Hz with 
a broadening of ±1.3 x 10^^^ x /o = ±0.05 Hz for the ring radius of 40 /im). For the first case, the fractional long 
term frequency fluctuation induced by the loss of one ion is 1.4 x 10^^^ x f{Mj) and is in the 10^^^ range for the 
larger ring. The choice of the Zeeman sublevel used for the clock operation results from a compromise between several 
shifts, including the Zeeman shift itself. In the next section we present these effects that are not directly related to 
the trapping electric field. 

IV. OTHER SYSTEMATIC SHIFTS AND THEIR DISPERSION 

A. Zeeman effect 

The selection rules for a quadrupole transition forbid AA/j — transition when the magnetic field lies along the 
direction of propagation of the clock laser. As the magnetic field is chosen along the trap symmetry axis to prevent 
the dispersion of the tensorial Stark shift, the use of AMj — transitions would require an angle between the trap 
axis and the direction of propagation of the laser. This option has to be rejected because of the large Doppler effect 
induced by the rf driven motion. As a second choice, the AM,/ = ±1 transitions are allowed when the laser propagates 
along the magnetic field direction. Indeed, the (S1/2, Mj = ±1/2 D5/21 Mj = Tl/2) transitions are only four times 
more sensitive to magnetic field fiuctuations than their {S1/2, Mj — ±1/2 T)^/2,Mj = ±1/2) counterparts. The 
summation of the frequency of the two transitions (S1/2, A^j — ±1/2 Mj = ±1/2) should cancel the first 

order Zeeman shift if the magnetic field is constant in time at every point of the ring. A field fluctuation SB induces 
an uncertainty and a line broadening on this transition of 2.2 MHz/G. From the frequency fluctuations evaluated in 
(36j for the JPL Hg+ microwave clock, and reached by a three layer magnetic shielding, one can infer magnetic held 
fluctuations of the order of 6 x 10"^ G, for a magnetic held of 0.05 G. 

The magnetic field has to be large enough to resolve the transitions between various Zeeman sublevels. In a Ca+ 
experiment, having a 1 kHz separation between neighboring transitions requires a magnetic field of 6 x 10^** G. If one 
assumes the relative stability demonstrated in the JPL microwave clock, the Zeeman shift fiuctuations are of the order 
of 0.013 Hz which result in a fractional uncertainty of 3 x 10"^^. If the same absolute stability can be reproduced, 
the Zeeman shift fluctuations reach 1.0 Hz which is equivalent to a fractional uncertainty of 2.5 x 10"^^. The same 
degree of frequency stability can be reached with less stringent conditions on the magnetic field stability using the 
hyperfine transition Mp = — > Mp = of the odd isotope ^'^Ca"'" [1^ [2^. Nevertheless, cooling and state detection 
with an odd isotope require more laser sources than with an even isotope and this complexity may not be acceptable. 

The Zeeman effect also reduces the stability of a single ion cloud. The main difference here is that the magnetic 
field has to be kept constant and homogeneous over the 40 fim (or 80 /im) of the ring diameter, compared to the 1 
fim scale associated to a single ion. This is still a favorable condition compared to the several centimeters long cloud 
used at JPL. 



B. Blackbody radiation shift 

The blackbody radiation shift (BBR shift) is the Stark effect induced by the thermal electric field radiated by the 
vacuum vessel and every part inside it. For calcium ions, very precise theoretical calculations [s^ for the polarizabilities 
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implied in the BBR shift aUow to know the shift to better than 3%: SfsBR = 0.38(1) Hz. This result assumes that 
the radiated field is isotropic and that the temperature of the emitting surface is 300 K. As the BBR shift scales 
with r^, thermal fluctuations can be detrimental to high precision clocks. A 10 K uncertainty keeps the frequency 
uncertainty and eventual broadening (less probabe as the time scales involved are very long) at the 0.05 Hz level, 
which is negligible in the context of our clock project. In conclusion, for temperature fluctuations smaller than 10 K 
a long term instability lower than 10^^^ can be reached. 

C. Quadrupole shift 

Another effect well known in single ion optical clocks is the quadrupole shift induced by the gradient of the local 
electric field coupled to the electric quadrupole of the Ds/^state (a S1/2 state has none). Thanks to a recent experiment 
[37} , in agreement with precise theoretical calculations ^sj, the quadrupole of the Ca+ D5/2 state is known to better 
than 1%: 0(D5/2)=1.83(1) ccq (atomic unit). The coupling with the rf electric field gradient gives rise to sidebands in 
the clock transition spectra that are well separated from the central band /o, so that the dc electric field gradient is 
the only relevant contribution for this shift. The dc trapping field has a quadrupole profile and using the notation of 
[33 |. its gradient is 2 A = — mw^/ {2q). For a magnetic field parallel to the trap symmetry axis and for uJz/2n — 1 MHz, 
the quadrupole shift is 

Sfg = 1.0 X (3il/^ - 35/4) Hz. (19) 

For the Zeeman sublevels Mj — ±1/2, the quadrupole shift induced by this dc trapping field is 8.0 Hz and does not 
depend on the position of the ions. If the axial trapping voltage is well controlled, this shift should not reduce the 
performance of the clock. It is not the case for the shift induced by extra dc fields which can build up in the trap, for 
example by neutral atom deposition. These fields can easily induce shifts of the order of 1 Hz from day to day operation 
and are well known from single ion clock operation. Two methods are used to reduce their day to day fluctuations: 
ionizing the neutral beam far from the clock operation area before shuttling them and/or using photoionization 39, 4^ 
which requires a far smaller flux of neutrals than the traditional electron bombardment method. If required, the shift 
itself can be compensated by the combination of three measurements on different Zeeman sublevels [4l| . The method 
based on frequency measurements with three orthogonal directions of the magnetic field, used for single ion standard 
BUI]; '^^^ ^ot be applied here because of the tensorial Stark shift dispersion fsee llll B"|) . All together, the uncertainty 
and fluctuations induced by this quadrupole shift must be made smaller than 0.04 Hz to keep the long term fractional 
frequency fluctuations smaller than 10"^^ . 

D. Misalignment 

If the propagation axis of the clock laser does not match the symmetry axis of the ring, the ions do not see the 
same laser power. The induced variations in the excitation probability contribute to the noise of the excitation 
detection and reduce the short term stability of the clock. Assuming a single Rabi pulse to interrogate the clock 
transition and a maximum excitation probability on resonance (^IlT = n with the Rabi frequency of the atom- laser 
coupling and T the pulse duration), the variation in the excitation probability for a detuning giving the maximum 
sensitivity is 6Pe{T) = 0.726^1^/^1^. If one considers a Gaussian laser beam of waist wl and intensity profile 
lL{r) = h exp(— 2r^/t«|), one can connect the variation S^l to a misalignment Sr. The dispersion in the probability 
of excitation is then 

ArSv 

sp^m ^o.m^^. (20) 

The quantum projection noise 1/2\/N j2l| is a fundamental limit which can only be beaten by entanglement and 
squeezing methods [l^. Making SPe{T) small compared to this noise is sufficient to make misalignment effects 
negligible. An easy solution is to make the laser waist a lot larger than the ring radius R. Already with wl — 2R, a 
geometrical precision of Sr < 0.27R is enough to keep SPe{T) < 0.1 and makes the contribution to the noise negligible. 

V. CONCLUSION 

The shift and uncertainty budget presented in table U shows that the rf is responsible for the major shifts and 
broadenings. The total broadening of the transition is smaller than 1 Hz, which is smaller than or comparable to 
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TABLE I: Uncertainty budget for the frequency transitfon of \Si/2,Mj = ±1/2) \D^I2,Mj = Tl/2) in '*°Ca+, 
based on a ring in an octupole linear trap with ijJz/Iti — 1 MHz, i7/27r=20 MHz, and a rf electric field such that 

R — 20/im (10 ions) or 40/im (20 ions), as given in the table. 
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the spectral broadening induced by the finite time excitation on the clock transition, depending on the Rabi pulse 
duration [2^. Therefore, the broadening due to the rf induced shifts does alter the short term stability and having 
10 or 20 ions instead of one effectively results in a gain in stability by \/N. 

Furthermore, the fluctuations of the rf induced effects over long time scales are not the limiting factors for the long 
term stability of the clock. Indeed, a single ion standard encounters the same limitations concerning uncertainty and 
long term stability, due to coupling to the local electric and magnetic fields. The challenge for a ring is to keep the 
stability constrains over larger spatial scales than for a single ion. 

As we have seen earlier, increasing the number of trapped ions allows to further reduce the long term fluctuations 
induced by ion loss, which scale like 5N/N^ and to lower the short term instability by increasing the signal to noise 
ratio. In this article, we limited our study to 10 and 20 ions ring to keep the temperature compatible with clock 
operation conditions while the cooling lasers are off. An alternative to a sequential operation that could allow to 
work with larger rings is to sympathetically cool the calcium ions by ions with a different mass (LC ions). Our MD 
simulations on 10 calcium ions show that conditions for axial and radial decoupling can be found, allowing the LC 
ions to reach their Doppler limit on the axial motion while calcium axial temperature fluctuates between 1 and 3 mK. 
These studies have to be extended [33*1 to larger samples to demonstrate a gain compared to the simple one species 
ring we have considered in this article. Such large cold rings could find an interest in metrology and/or quantum 
information or simulation. With today's state of the art in magnetic field stabilisation, the ion loss is not the limiting 
effect for the frequency stability of a clock based on a ring of ions and this configuration can offer the possibility to 
test many particles system for metrology. 
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